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Abstract 

Based on extensions of the grand-canonical Quantum Monte-Carlo algorithm to 
incorporate magnetic fields, we provide numerical data confirming the existence of 
a Kosterlitz-Thouless transition in the attractive Hubbard model. Here, we calcu- 
late the temperature derivative of the superfiuid density, ^^g^^^^ ? to pin down the 
transition. The latter quantity is obtained from the difference in internal energy 
of systems which differ by a phase twist 7r/2 in the boundary condition along one 
lattice direction. Away from half-band filling, ^^^"J^^^ shows a response which in- 
creases with lattice size at the transition temperature. In contrast, such a signal is 
not observed for the case of a half-band filling. 
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Numerical simulations of the attractive Hubbard Hamiltonian indicate that 

this fermionic lattice model has an s-wave superconducting ground state in two- dimensions. 
A finite-size scaling analysis (up to one hundred sites) of the s-wave pair-pair correlation 
functions has provided some evidence that the model has a finite-temperature Kosterlitz- 
Thouless (KT) [Q transition to a superconducting state away from half-band filling 0. 
Based on extensions of the grand-canonical Quantum Monte-Carlo (QMC) algorithm, we 
present numerical data which clearly confirm the conjecture of a KT transition. Our 
approach relies on the calculation of the temperature derivative of the superfluid density, 
^^aff^^^ , where (3 is the inverse temperature. In the framework of such a KT transition, 
Ds{P) shows a universal jump at the transition temperature Hence, ^^g^^^'* in the 
thermodynamic limit behaves like a Dirac delta-function at the transition temperature, 
and even on finite-sized clusters is expected to yield a very clear signal. Compared to 
previous estimations of the T = superfiuid density , via fiux quantization P or linear 
response [0, the method described here leads to numerical results which exhibit a peak 
at TxT which grows as the size of the system increases. This provides a clear signal of 
the KT transition. Such signatures of the KT transition have already been observed in 
classical ^ as well as in quantum XY models 0. In this letter it is shown that, for the 
attractive Hubbard model, ^^^"J'^^ is given by the difference in internal energy of systems 
which differ by a phase twist 7r/2 in the boundary condition along one lattice direction. 
The internal energies are extracted from finite temperature quantum Monte-Carlo (QMC) 
data on cluster sizes ranging from 4 x 4 to 8 x 8. Comparison of the data for different 
lattice sizes, strongly supports the existence of a KT transition away from half band filling. 

The attractive Hubbard hamiltonian we consider is given by: 



tyici c^_^. + 
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+ f^El^^T-^) -/^EKt + ^u). u <q. (1) 



Here, cl^ creates an electron with z-component of spin a on lattice site Ex^y are the 
lattice vectors of unit length, $o = — is the fiux quantum, and L is the linear length of 



the square lattice. The fermionic operators obey periodic boundary conditions in both 
lattice directions. The above Hamiltonian thus describes an attractive Hubbard model 
on a torus. The phase factors in the kinetic energy term account for the presence of a 
magnetic field of flux $ threaded through the center of the torus. Through a canonical 
transformation, one may remove the phase factors in (|I]). However, under this canonical 



transformation, the x-boundary acquires a phase twist exp (27ri$/$o) 1131 • Hence, the 
threaded magnetic field acts merely as boundary effect. 
Let us now consider the free energy, 

= -^InTre-'^^^*), (2) 

P 

and expand it up to second order in $/$o around $ = 0. This yields: 



= F{^ = 0) + {—j D,iP) + OU — ) ) where 



(^<K,>+l^dr<J,{T)J,{0)>y 

Jx = —ity^ (cl Ci'.ff „ — cl ^ C7„) 

= -ty(cl G,. ^ + c^J. (3) 

Here, Ds{i3)/2tt'^ corresponds to the superfiuid density in units where the hopping t = 1 



rT| . The above equation may be transformed to give: 

|(W4)-W4 = 0)) = (^)^|(M+o((|)'), (4) 

The left hand side of the above equation is just the difference in internal energy: £'($) — 
= 0). The right hand side contains the term we are interested in; namely the 
temperature derivative of the superfiuid density. For an infinite lattice, the jump in the 
superfiuid density at T^t implies that 

— — (5 (T - Tkt) ■ (5) 

For a finite-size lattice, the delta-function peak scales as the number of sites P]. 

To determine which values of $ are most appropriate for our calculations, we consider 
the T = limit. At zero temperature, the attractive Hubbard model is expected to 



show off-diagonal long-range order (ODLRO), and thus ffux quantization |T^. Hence, the 
quantity 

Eo($) - Eoi^ = 0) = lim lim (F($) - F($ = 0)) (6) 

scales in the thermodynamic limit to a periodic function of period $o/2. Furthermore, a 
non- vanishing energy barrier is to be seen between the flux minima at $ = and $ = $o/2 
[llO| , Thus, the maximal response in Eo{^) - Eo{^ = 0) should lie at $ = $o/4. The 
superfluid density, which corresponds to the curvature of the -E'o('^') — Eq{^ = 0) function 
at $ = or $ = $o/2 [0, 0], may be approximated by: 

^'^^r ~ ^o('^ = "^0/4) - ^o(<f = 0) = ^o($ = $o/4) - ^o($ = $o/2). (7) 

The above relations hold for the bose condensation of free charge two bosons on a square 
lattice. Numerical evidence for the validity of equation for the quarter-filled and 
half-filled two-dimensional attractive Hubbard model may be found in reference 0. 
Due to the above considerations, we have computed the quantity: 

AE{R) = = $o/4) - = R) where R = R = %/2 (8) 

as a function of temperature and lattice size. AE{R) yields at finite-temperature infor- 
mation on ^^g^*-^'' d^) and at T = converges for increasing lattice size to the superfluid 
density (0). Considering both R = i? = $o/2 yields a cross-check for finite size ef- 
fects since both values of R converge to the same results only in the thermodynamic 
limit. As mentioned previously, the internal energies E{^) were calculated through the 



use of the finite temperature QMC algorithm [1^, We carry out a discrete Hubbard- 
Stratonovitch (HS) transformation of the on-site interaction term in (||). After integration 
of the fermionic degrees of freedom, one obtains for the attractive Hubbard model: 

Z = Tre-'^^(*) = 5^(detM(s))' (9) 

s 

where s denotes a HS configuration and M is an L x L matrix. For $ = Oor$ = $o/2 
(i.e. antiperiodic boundary in the x-direction), M is real and thus no sign problem 
occurs. On the other hand, for $ = $o/4, M is complex and thus, (detM(s))^ is equally 
complex. Hence, this value of the flux yields a sign (or in our case a phase) problem. 



This phase problem inhibits us of reaching very low temperatures. Another problem 
we encountered, is that of fixing the chemical potential. Since we are interested in the 
difference of two internal energies at a given particle number (^), the chemical potential 
has to be determined very precisely. Apart from the half-filled band case where /i = a 
brute force search for the desired chemical potential is numerically too expensive. We thus 
carried out at least three simulations with chemical potential corresponding approximately 
to the desired filling and fit the energy data to the form: A + Bx + Cx"^. The energy 
corresponding to the desired filling is then extrapolated from the fit. Thus, one data point 
in the AE{R) curve corresponds to at least six independent simulations. 

Figure 1 plots AE{R) for both considered values of R (see equation (^)), and for the 
half-filled attractive Hubbard model at U/t = —4 as a function of the temperature. The 
data points at T = were obtained with the projector QMC algorithm (PQMC) p!5| . 
Here, we consider three lattice sizes: 4x4, 6x6 and 8x8. At exactly half-band filling, the 
attractive Hubbard model may be mapped onto the repulsive Hubbard model by carrying 
out a particle-hole transformation in one spin sector. Under this transformation, one 
notices that the s-wave pair-pair (charge-charge) correlations of the attractive model map 
onto the transverse (z-component) spin-spin correlations of the repulsive model. In the 
strong-coupling limit, the repulsive Hubbard model may be mapped onto the Heisenberg 
model which shows no KT transition [|TB|. Hence, the strong-coupling attractive Hubbard 
model at half-band filling shows no KT behavior. In the intermediate coupling regime, 
the numerical data of figure 1, equally display no sign of a KT transition. Consider the 
R = $o/2 curves. For this value of R, there is a bump appearing approximately at 
T ~ 0.2t. However, this structure fails to scale with increasing lattice size but rather 
appears to be suppressed. A similar, but less pronounced structure may be seen in the 
two largest lattice sizes of the R = curves. 

Let us now turn to three-quarters band filling. Our QMC data are plotted for R = 
{R = 0.5) in figure 2a (2b). For the considered coupling {U/t = —4) and band-filling 
(< n >= 0.75), we were unable to reach inverse temperatures higher than jSt = 7 due to 
severe phase problems. Again, the data points at T = were obtained with the PQMC 
algorithm |T^. Figures 2a and 2b show a striking difference between the 4x4 and 8x8 
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clusters. The slope of the AE{R) curves become very large for the 8x8 cluster size in 
the temperature region 7 < pt < 5 ^7^. Finally, the available data, are compatible with 
the onset of a Dirac delta- function type response in AE{R) for both considered values 
of R at approximately T ~ O.lt. This temperature value compares favorably with the 
KT transition temperature obtained from a finite-size scaling analysis of the pair-pair 
correlation functions 0. 

In conclusion, extensions of the grand-canonical QMC algorithm have allowed us to 
calculate the temperature derivative of the superfiuid density for the half and three- 
quarters filled attractive Hubbard model at U/t = —4. Comparing the numerical results 
for both band fillings provides convincing evidence that the three-quarters filled attractive 
Hubbard model at U/t = —4 has a KT transition to a superconducting state. In contrast, 
no signature of the KT transition is found at half-band filling. 
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Assaad would like to thank the DFG for financial support under the grant number Ha 
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Figure captions 

Fig. 1 = $0/4) — = i?) as a function of temperature for the half-filled 
attractive Hubbard model at U/t— —4. Here, we consider three lattice sizes: 
4x4, 6x6 and 8x8. The three upper (lower) curves correspond to R — 
(R — O). The solid hues are a sphne fits to the data. The data points at T = 
were obtained with the PQMC algorithm. 

Fig. 2 a) — $0/4) — = 0) as a function of temperature for the three- 
quarters filled attractive Hubbard model at U/t = —4. The data points at 
T = were obtained with the PQMC algorithm. The solid lines are spline fits 
to the data. 

b) Same as figure a) but we consider: = $o/4) - E{^ = ^o/2). 
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